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We investigate the collective behavior of orbital angular momentum in the spin ferromagnetic state 
of a Mott insulator with t^ g orbital degeneracy. The frustrated nature of the interactions leads to an 
infinite degeneracy of classical states. Quantum effects select four distinct orbital orderings. Two 
of them have a quadrupolar order, while the other states show in addition weak orbital magnetism. 
Specific predictions are made for neutron scattering experiments which might help to identify the 
orbital order in YTiC>3 and to detect the elementary orbital excitations. 
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Recently, "orbital physics" has become an important 
topic in the physics of transition metal oxides 0,0 . Or- 
bital degeneracy of low energy states and the extreme 
sensitivity of the magnetic bonds to the spatial orienta- 
tion of orbitals lead to a variety of competing phases that 
are tunable by a moderate external fields. A well known 
example is "colossal magnetoresistivity" of manganites. 
A "standard" view on orbital physics is based on pio- 
neering works developed and summarized in Ref. 
|^| . In this picture, long-range coherence of orbitals sets 
in below a cooperative orbital/ Jahn- Teller (JT) transi- 
tion; spin interactions on every bond are then fixed by 
the Goodenough-Kanamori rules. Implicit in this picture 
is that the orbital splittings are large enough so that 
one can consider orbital populations as classical num- 
bers. Such a classical treatment of orbitals is well jus- 
tified when orbital order is driven by strong cooperative 
lattice distortion. 

However, interesting quantum effects are expected to 
come into play when orbital physics is dominated by elec- 
tronic correlations. Mott insulators with threefold ti g 
degeneracy (titanates, vanadates) are of particular inter- 
est in this respect, as the relative weakness of the JT 
coupling for t2 g orbitals and their large degeneracy en- 
hance quantum effects. A picture of strongly fluctuating 
orbitals in titanate LaTiOa has recently been proposed 
on both experimental [[| and theoretical grounds . 

Titanates are very interesting: they show a nearly 
continuous transition from an antiferromagnetic state in 
LaTi03 with unusually small moment to a saturated fer- 
romagnetic one in YTi03 [||. Yet another puzzle is the 
almost gapless spinwave spectrum of cubic symmetry re- 
cently observed in YTi03 by Ulrich et al. |[o]]. This 
finding is striking, because the orbital state predicted 
by band structure calculations |ll| is expected to break 
cubic symmetry of the spin couplings. In principle, one 
may obtain isotropic spin exchange couplings by using 
this orbital state |jl2); however, this requires fine tuning 
of parameters, and rather small changes of the orbital 
state may reverse even the sign of the spin couplings [Jl0[ . 



To gain more insight into the origin of the different be- 
havior of orbitals in YTi03 and LaTiOs, we ask in this 
Letter the following question: What kind of orbital state 
would optimize the superexchange (SE) energy if we fix 
the spin state to be ferromagnetic as in YTi03? We find 
a surprisingly rich physics: (i) there are several "best" or- 
derings, some of them are highly unusual showing static 
orbital magnetism; (ii) the resulting spin interactions re- 
spect cubic symmetry in accord with experimental data; 
(iii) strong orbital fluctuations lead to anomalously large 
order parameter reduction. These states have been over- 
looked previously, as their stabilization is an effect of 
quantum origin. We also discuss a mechanism which sta- 
bilizes ferromagnetic state in YTi03. 

The t-2 g superexchange can, in general, be written as 

H SE = JseE W ) [{^ ■ S 3 + |) + where 
Jse = 4:t 2 /U represents the overall energy scale. The or- 
bital operators and K^' depend on bond directions 
7(= a, 6, c). Their explicit form in a d fog) system like 
the titanates is given in Ref. ||. If spins are all aligned 
ferromagnetically (i.e., Si ■ Sj = 1/4), the above equation 
reduces to the orbital Hamiltonian: 

ff-rb = j£4 7) - (1) 

m 

Hereafter, we use r\JsE as a unit of energy and drop 
a constant energy shift (= —tiJse)- The coefficient 
7*i = 1/(1—377) originates from Hund's splitting of the ex- 
cited t\ g multiplet via 77 = Jh/U. The orbital operators 

AlH' can be represented via a triplet t, = (a, 6, c)$ of con- 
strained particles (n ia + n ib + n ic = 1) corresponding to 
tig levels of yz, xz, xy symmetry, respectively. Namely, 

A[f = n. la n ja + n lb n jb + nfiji^i, + Wn.a'Ji,. (2) 

for the pair along the c axis. Similar expressions are ob- 
tained for the exchange bonds along the a and b axes 
by using in Eq. (Q) (6, c) and (c, a) orbiton doublets, re- 
spectively. The Hamiltonian ([!]) is the main focus of 
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this Letter; a theory of ferromagnetic magnons including 
spin-orbital effects will be presented in Ref. [^3| . 

It is useful to look at the structure of H or i, from differ- 
ent points of view (see also Refs. |7|,p|Jl4|] ) : (i) On a given 
bond, the operator Ay acts within a particular doublet 
of orbitals. Spinlike physics, that is, the formation of or- 
bital singlets, is therefore possible. This is the origin of 
orbital fluctuations in titanates. (ii) On the other hand, 
interactions on different bonds are competing: they in- 
volve different doublets, thus frustrating each other. This 
brings about a Potts model like frustrations, from which 
the high degeneracy of classical orbital configurations fol- 
lows, (iii) Finally, using the angular momentum opera- 
tors l x = i(c^b — wc) , l y = i(alc— c'a) , l z — i(b<a — aJb) 
of t2 g level, one may represent the interaction as follows: 

~\~(lxly)i(lylx)j ~t~ (lylx)i{lxly) j ■ (3) 

(For b and a bonds, the interactions are obtained by per- 
mutations oi l z ,ly,l x .) We observe a pseudospin I = I in- 
teraction of pure biquadratic form. Would / be a classical 
vector, it could change its sign at any site independently. 
A local "Z2 symmetry" and the associated degeneracy 
of the classical states tells us that angular momentum 
ordering, if any, is of pure quantum origin. 

The above points (i) — (iii) govern the underlying 
physics of the orbital Hamiltonian. By inspection of the 
global structure of Ay , one observes that the first two 
terms are definitely positive. However, the last two terms 
(which change the "color" of orbitals) can be made nega- 
tive on all the bonds simultaneously, if (i) on every bond, 
two particular components of U and lj are antiparallel, 
and (ii) remaining third components are parallel. For c 
bonds the rule reads as: U z lj Z and U x ljx are both nega- 
tive, while ly components are parallel. (In terms of or- 
bitons: Ci and Cj are in antiphase, ai and aj as well; but bi 
and bj have the same phase.) We find only two topolog- 
ically different arrangements [called (a) and (b)], which 




FIG. 1. Arrangement of the local quantization axes in 
states (a) and (b). Each state has four sublattices denoted by 
numbers. Arrows show a direction of the quantization axes 
at each site. They represent also a snapshot of local correla- 
tions of orbital angular momentum: on every bond, two out 
of three components l a are correlated antiferromagnetically. 



can accommodate this curious mixture of "2/3 antiferro" 
plus "1/3 ferro" correlations (see Fig.l). 

For technical reasons, and also to simplify a physical 
picture, it is useful to introduce new quantization axes. 
This is done in two steps. First, we introduce local, sub- 
lattice specified quantization axes (see Fig.l): 

1 : (x, y, z) -> (x, y, z), 2 : (x, y, z) -> (-x, -y, z), 
3 : (x, y, z) -> (-x, y, -z), 4 : (x, y, z) -> (x, -y, -z). 

After corresponding sign transformations of li a and or- 
bitons, one obtains a negative sign in front of the last 
two terms of Ay\ thereby "converting" the interactions 
in a new axes to a fully ferromagnetic one. From now 
on, a sublattice structure will not enter in the excitation 
spectrum. From the above observations, it is also clear 
that all the components of I are equally needed to opti- 
mize all three directions. We anticipate, therefore, that 
the cubic diagonals are "easy" (or "hard" ) axis for I fluc- 
tuations/orderings (recall that the Hamiltonian has no 
rotational symmetry). The second step is then to rotate 
quantization axes so that the new z axis (z) is along one 
of the cubic diagonals (say, [111])- This is done as fol- 
lows: U = R U (and t = R t for orbitons). The matrix 
R = R z (n/4)Ry(6 )R z (-ir/4) with tan(9 = V2. R a {9) 
describes the rotation around a axis by angle 9. 

By the above transformations, the orbital Hamiltonian 
obtains well structured form and can be divided into two 



parts as H, 
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H 



I sing ' 



Here, Hj sinq represents 



a "longitudinal" part of the interaction: 
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Htrans is responsible for fluctuations of transverse com- 
ponents l x , ly and quadrupole moments of various sym- 
metry. Its explicit, somewhat lengthy form will be pre- 
sented in Ref. fll3|] . In terms of a, b, c orbitons, Hi S i ng 
promotes a condensation of an appropriate orbiton, while 
Htrans gives a dispersion of orbitons and their interac- 
tions. 

Order and excitations. — Now we are ready to discuss 
possible orbital orderings. Equation (Q) makes it explicit 
that there are two deep classical minima corresponding to 
two different "ferro-type" states: (I) quadrupole ordering 
(classically, k z — 0) and (II) magnetic ordering (k z — 1, 
classically). We notice that the last, magnetic term in 
Eq. (^) is generated by quantum commutation rules when 
we rotate H rb\ this makes explicit that the Z2 symmetry 
is only a classical one and emphasizes the quantum origin 
of orbital magnetism. We see below that states I and II 
are degenerate even on a quantum level. Noticing that an 
arbitrary cubic diagonal could be taken as z and having in 
mind also two structures in Fig.l, one obtains a multitude 
of degenerate states. This makes, in fact, all the orderings 
very fragile. 
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To be more specific about the last statement, we 
calculate the excitation spectrum. Introduce first a 
quadrupole order parameter Q — ((l 2 + ly)/2 — I 2 ) = 
(rig— (n.a+nt)/2), and the magnetic moment mi — (l z ) — 
i(b'a — d^b) (we set [is = 1)- State I corresponds to the 
condensation of the orbital c (an equal mixture of the 
original a,b,c states), while the ordering orbital in the 
state II is an imaginary one (a — ib) /y/2. Let us call the 
remaining two orbitals by a and /3, which are the exci- 
tations of the model. Specifically, a = a, = b in state 
I, and a = [a + ib + (1 + i)c]/2, = [b - ia - (1 - i)c]/2 
in state II. Note that latter doublets can be regarded as 
magnons of orbital origin, describing fluctuations of the 
orbital magnetism. Remarkably, in terms of these dou- 
blet excitations, the linearized Hamiltonian in a momen- 
tum space has the same form for both phases. Namely, 

H ow = Efc [n a k + n /3k + { l(7i + 72K"_fe 

+|(7i - 72)/%/?_g - 73«j/3_g + H.c.}} , (5) 

where 71 = (c x + c y + c z )/3, 72 = \fi(c y - c x )/6, 73 = 
(2c z — c x — c y )/6, and c a = cosfc Q . These expressions 
refer to state I. For state II, one should just interchange 
72 and 73. However this does not affect the spectrum of 
the elementary excitations, which is given by co±(k) = 
s/l - (71 ± k) 2 , where k = y/i'£ +73. 

Using matrix Green's functions for the (a, (3) doublet, 
we may calculate various expectation values: (i) quantum 
energy gain: Eq = —0.214 for all the phases; (ii) num- 
ber of bosons not in the condensate: (n a + np) = 0.54. 
Thus, the order parameters are strongly reduced by fluc- 
tuations. We obtain Q = 0.19 (mi = 0, of course) in state 
I, and mi = 0.19 accompanied with small quadrupole mo- 
ment Q — —0.095 in state II. To visualize the orbital pat- 
terns, we show in Fig. 2 the electron density calculated 
including quantum fluctuations. 

The anomalous reduction of order parameters is due to 
the highly frustrated nature of the interactions in Eq. ([!]). 
A special, nonspinlike feature of all orbital models is that 
orbitals are bond selective, resulting in a pathological de- 
generacy of classical states. This leads to soft modes [ob- 
serve that oj±(k) is just flat along (0, 0, tt) and equivalent 
directions]; such soft modes were first noticed in Ref. 
fL5| . Linear, "spin- wave-like" modes are, however, not 
true eigenstates even in the small momentum limit, as 
orbital models have no continuous symmetry. Therefore, 
an orbital gap must appear (see Ref. Jl6| and, in par- 
ticular, the discussion of the "cubic" model in Ref. M). 
In the present model, we have estimated the orbital gap 
A = 0.57 using a single mode approximation ]l7| . 

Angular momentum fluctuations can directly be 
probed by neutron scattering measurements. We have 
calculated the structure factor S(q,u>) — -Imx(<f, oj), 
where xiQ^) is the orbital angular momentum suscep- 
tibility. We first obtain x(q, oj) in a rotated basis. A 



noncollinear, four sublattice orbital order generates in 
S(q, oj) also an additional terms x(q + <fi, oj) with shifted 
momenta. For the arrangement (a) q\ — (n, 0, tt), q? = 
(-7T, 7T, 0), <?3 = (0,7r,7r), while for the structure (b), one 
has <fi = (tt , 7r , 7r) , g*2 = (tt, 7r, 0), <f 3 = (0, 0, tt). The re- 
sults are shown in the right panels in Fig. 2. Neutron 
spectroscopy at higher than spin-wave energies are nec- 
essary to detect "orbital magnons." If one of the states 
11(a), (b) is realized, a static Bragg peak of orbital ori- 
gin must show up. Reexamination of form-factor data 
by Ichikawa et al. using t2 g spin densities shown in 
Fig. 2 is also desirable. 

A remarkable feature of all the above phases is their 
robust property to provide exactly the same spin cou- 
plings in all cubic directions, hence resolving the puz- 
zling observation by Ulrich et al. jl0|]. A small spin- 
wave gap A s in YTiO*3 further supports the theory. In 
state I (a), we obtained a spin-orbit coupling A- induced 
gap A s ~ A 2 /J S e, where A ~ X 2 /3riJ S E- With 
A ~ 15 — 20 meV and Jse ~ 40 — 50 meV, one obtains A s 
below 0.1 meV consistent with experiment. The gap is 
small due to high symmetry of the orbital ordering. On 
the other hand, a large classical gap A s ~ A is obtained 
in states I (b), II (a), (b). Thus, the data of Ref. 0] 
supports state I (a) in YTiC>3. 

We discuss now our findings in a context of the 
interplay between different spin states in titanates. 
It is noticed that the quantum energy gain Ep = 
—0.21AriJsE obtained above is only slightly less than 
the energy of a spin AF/orbital disordered phase, Eaf — 
—0.33 \^ T\J se Jl9[ stemming from a composite spin- 
orbital resonance. We argue now that the distortion of 
Ti-O-Ti bond angle due to a small size of Y ion further 
reduces the difference AE = Eaf — Ep and stabilizes 
the ferromagnetism of YTiC>3. This distortion (i) reduces 
hopping t = t cos 9 ||; hence, Jse = J S ° E cos 2 9, and (ii) 
creates finite transfer t' = t' sin 9 between t^g and e g or- 
bitals [see inset of Fig. 3 (a)]. Here t' Q /t = td pa /td im ~ 2, 
and idpo-(ir) is the <j(tt) bond transfer between Ti 3d and O 
2p orbitals. The channel (ii) generates unfrustrated fer- 
romagnetic interaction (caused by Hund's coupling be- 
tween t2 g and e g electrons in the excited state) with 
J' = -|4°i(L//C/) 2 7 ? sin 2 9. Here U = U + A cr with 
A cr being a cubic crystal field splitting. As a result, we 
obtain AE = [-0.11 cos 2 9 + 2{U/U) 2 ij sin 2 9}J ( ° E . As 
shown in Fig. 3 (a) , AE reduces with decreasing 9 and a 
transition from AF to ferrostate occurs at critical angle 
136°. A ferromagnetic state with orbital order 



n(0) 



is expected to be further supported by small JT energy 
gain. By introducing 5Ejt = — 0.04J 



(0) 



SE 



~> —3 meV, we 
obtain 9 C = 146° as observed in titanates We notice 
that concominant JT distortions associated with such a 
small SEjp are expected to be very weak. 

Finally, it is stressed that orbital orderings are found 
to be weak; hence, large scale orbital fluctuations must 
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be present. This implies strong modulations of the spin 
couplings, both in amplitude and sign, suggesting a pic- 
ture of "fluctuating exchange bonds", where the mag- 
netic transition temperature reflects only a time average 
of the spin couplings. The intrinsic quantum behavior of 
t 2g orbitals in titanates may lead to a scenario, outlined 
in Fig. 3 (b), for an explanation of the puzzling "smooth" 
connection between ferromagnetism and its antipode, the 
Neel state. This proposal can be tested experimentally: 
we expect almost continuous magnetic transition under 
pressure or magnetic field that conserves the cubic sym- 
metry of spin exchange constants. 

To conclude, the ferromagnetic phase of the t^g su- 
perexchange system has several competing orbital order- 
ings with distinct physical properties. The elementary 
excitations of these states arc identified. We find that 
Ti-O-Ti bond distortion favors the ferromagnetic state. 
Theory well explains recent neutron scattering results in 
YTiC>3. An orbital driven quantum phase transition in 
titanates remains a major challenge for future study. 
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FIG. 2. t2 g electron density (left) and intensity of the or- 
bital contribution to the structure factor S(q,u>) (right). Up- 
per panel: Quadrupole ordered state I (a). Because of the 
more spatial anisotropy of electron distribution, this state 
may be preferred by electron-lattice coupling. "Hot spot" 
at (tt, n, 7r) point is due to coherency factors stemming from 
four sublattice structure. Lower panel: Magnetic state II (a). 
Shown in the inset are the expected orbital Bragg peak posi- 
tions. 
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FIG. 3. (a) Energy difference (in units of Jg'g) between AF 
and ferromagnetic states as a function of 9. The solid line is 
obtained within the SE model. Parameters are r\ = 0.12 and 
U/A cr = 2.5. The broken line includes small JT energy gain 
in the orbitally ordered ferromagnetic state, (b) Qualitative 
picture of the evolution of interactions in the spin channel. 
An average, static component of the spin exchange "constant" 
seen by coherent spin-waves is represented by a solid curve. 
Its sign depends on local correlations of orbitals. To the right 
of the critical point (QCP) these correlations are more an- 
tiferromagnetic, supported by noncollinear orbital orderings. 
To the left, the genuine ground state of the ti g system on 
cubic lattice, an orbital disordered state supporting spin AF, 
is stabilized. In the proximity area, a fluctuating part of the 
overall SE interaction dominates, and separation of the spin 
and orbital degrees of freedom might no longer be possible. 
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